Identifying and Indexing Icosahedral Quasicrystals from Powder Diffraction Patterns 
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We present a scheme to identify quasicrystals based on powder diffraction data and to provide a 
standardized indexing. We apply our scheme to a large catalog of powder diffraction patterns, in- 
cluding natural minerals, to look for new quasicrystals. Based on our tests, we have found promising 
candidates worthy of further exploration. 

PACS numbers: 61.10.Nz, 61.44 Br 



Quasicrystals are solids whose lattices exhibit a rota- 
tional symmetry, such as five-fold symmetry, that is for- 
bidden for periodic crystals. Because the atoms are 
arranged quasiperiodically, the structure's symmetry is 
not subject to the usual mathematical restrictions of 
crystallography. Quasiperiodic translational order has 
physical consequences. For example, since electrons and 
phonons in quasicrystals do not encounter a periodic 
potential, quasicrystals have unusual resistive and elas- 
tic properties, and these have been exploited in several 
applications. M To date, known quasicrystals have been 
found by serendipity or by probing stoichiometric vari- 
ations around other already-known quasicrystals. Fur- 
thermore, all known quasicrystals are synthetic; no nat- 
ural quasicrystal has ever been identified. A more sys- 
tematic way to search for quasicrystals, including natu- 
ral quasicrystals, is desirable, and one way is to search 
through collections of diffraction data. Although a two- 
dimensional electron diffraction pattern would immedi- 
ately show a quasicrystal's salient forbidden symmetry, 
no large collection of such patterns exists. However, 
a collection of over eighty thousand powder diffraction 
patterns in digital form, the Powder Diffraction File 
(1CDD-PDF), is published by the International Center 
for Diffraction Data. The catalog contains synthetic in- 
organic and organic phases, as well as nine-thousand min- 
eral patterns. Because the powder diffraction pattern of 
a material averages over all orientations, only the mag- 
nitude (and not the direction) of the scattering vector in 
reciprocal space is preserved, and a quasicrystal's distinc- 
tive non-crystallographic symmetry cannot be observed 
directly. A priori, it is unclear if quasicrystals can be 
identified from their powder diffraction patterns alone. 

In this paper, we present a method to identify, classify 
and index icosahedral quasicrystals based solely on their 
powder diffraction patterns. We apply the method to 
the ICDD-PDF and hnd the best-fit quasicrystal candi- 
dates. Our automated procedure picks out the known 
quasicrystals in the ICDD-PDF and, for the cases in 
which indexing has been published, produces the same 
indices. Although the remainder of the catalog is sup- 
posed to consist of periodic crystals, it was assembled 



over many years, including decades prior to the 1984 dis- 
covery of quasicrystals. Consequently, it is conceivable 
that the catalog includes some quasicrystals that were 
never identified as such. Based on our studies, we report 
promising materials. 

The diffraction pattern of an ideal three-dimensional 
quasicrystal consists of Bragg peaks located on a lattice 
given by Q — 5Zi=i n J>ii where the bi are basis vectors 
pointing to the vertices along the six five-fold symme- 
try axes of a regular icosahedron in three dimensions, 
and the are integers that index each vector. The 
quasicrystal formula is similar to that for a crystal ex- 
cept that the number of basis vectors is greater than the 
number of dimensions (three), a consequence of the non- 
crystallographic symmetry. We choose the bi following 
the convention in jj| : 

& 1 = (l,r,0), b 2 = (r,0,l), 6 3 = (0,l,r), (1) 
& 4 = (-1,t,0), b 5 = (r,0,-l), 6 6 = (0,-Lr) 

where r is the golden ratio, (1 + v5)/2- An equiva- 
lent way to index the position of the three-dimensional 
reciprocal-space scattering vector is to use a scheme anal- 
ogous to crystallographic Miller (h k I) indices. Six inte- 
ger indices grouped into pairs describe the distance along 
Cartesian basis vectors, Q = (h + h'r)x + (k + k'r)y + 
(I + Vt)z. The (h/h' k/k' l/V) indices are permutations 
of the ny. 
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k' 



m 



■ «4, 



hi = 712 + m 

I = n 2 - n 5 , 



n 3 
n 3 



(2) 



Because the {h/h' k/k' l/V) indices express distances 
along Cartesian axes, the advantages of orthogonal co- 
ordinate axes are conveniently recovered. Also, associ- 
ated with every Q is a vector Q±, constructed from an- 
other integer linear combination of the same basis vec- 
tors: Q± = (h'-hT)x+(k'-kr)y+(l'-lT)z = £Li nSi, 
where 

bi = (-r,l,0), bi - (1,0, -t), bi = (0, -r, 1), (3) 
6f = (r,l,0), bi = (1,0, r), 6^ = (0,r,l) 
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In powder diffraction patterns, the reciprocal space is col- 
lapsed to one dimension, where all vectors with the same 
magnitude \Q\ are degenerate. The magnitude \Q\ = Q 
is the reciprocal of the d-spacing, the quantity listed 
in the entries of the ICDD-PDF. Q 2 and \Q ± \ 2 = Q\ 
can now be expressed as integer linear combinations of 
only two basis vectors, whose lengths are related by the 
golden ratio: Q 2 oc N + tM and Q\ cx Nt — M, where 
N = 2£)f =1 n 2 and M = h' 2 + k' 2 + l' 2 +2(hti + kk' + W). 
Each peak in a quasicrystal powder diffraction pattern 
can be indexed by the two integers N and M. The entire 
diffraction pattern can be scaled by a factor of r 3 along 
the Q direction, and the support will remain the same, 
a consequence of the self-similarity of the pattern. That 
is, for each lattice vector at position Q 2 = (N,M), there 
is another lattice vector of similar intensity at position 



(N',M') 



Using the relation r = r + 1 , 



(N', M') =N' + M't = t 6 (N + Mr) = (4) 
(5M + 8N) + (8M + 13JV)t = (5M + 8N, 8M + 13N) 

Although three other different indexing schemes persist 
in the literature to describe icosahedral quasicrystal pow- 
der patterns, [0, [| [| they are all analytically equivalent 
to the convention given here. 

There are three distinct icosahedral reciprocal lattices: 
simple icosahedral (SI), face-centered icosahedral (FCI) 
and body-centered icosahedral (BCI). The SI, FCI and 
BCI reciprocal lattices correspond to real-space lattices 
obtained by, for example, placing identical "atoms" at 
each lattice point of a simple-hypercubic, face-centered 
hypercubic and body-centered hypercubic lattice in six 
dimensions, respectively, and projecting down to three 
dimensions. (In |Q, the SI, FCI and BCI lattices are 
referred to as P*, I* and F*, respectively.) For these 
primitive lattices, the intensity is given by:|| 



sin(Q ± /2) 



Q±/2 



(5) 



Note that I increases as Q±^ — ► so that bright peaks — 
the ones likely to be observed experimentally — have small 
Q±. Eq. (||) does not account for chemical and other ef- 
fects that modulate the intensities in the diffraction pat- 
tern of a real material, but it provides a guide as to which 
peaks should be observed when testing real patterns. 

Testing patterns. The first step in testing a real powder 
pattern is to find the best possible peak-by-peak match 
between that pattern {Qi} and the perfect quasicrystal 
template pattern Then, various statistical tests can 
be applied to compare that match to what is expected for 
a true quasicrystal. 

If we were to match to a periodic crystal template pat- 
tern, the test would be more straightforward. For any 
finite interval of Q over which a pattern is measured, a 
real crystal has a finite number of diffraction peaks. One 



also has good a priori estimates for the nonzero mini- 
mum distance between nearest peaks, which is set by the 
lattice constant. These features, useful in matching the 
template to the real pattern, allow a unique indexing. 

For the quasicrystal, the process is more complicated 
because the number of peaks in any finite interval of the 
ideal pattern is infinite (the pattern is dense) and, conse- 
quently, there is no uniquely defined lattice constant or 
indexing. How does one sensibly match a real pattern, 
with a finite number of peaks in a given interval AQ, 
to a quasicrystal template pattern, with a dense set of 
peaks in that same interval? One must account for the 
intensity of the peaks, not just their Q. Although the 
quasicrystal powder pattern is dense, most of its peaks 
within any finite interval of Q have a large Q± and are 
predicted by Eq. (Q) to be too dim to distinguish from the 
noise present in any real experiment. So, instead of find- 
ing the template which comes closest to the observed 
Qi (i.e. minimizing \Qf — qf,\ alone), we instead mini- 
mize \Q\ — qf,\/Ii>, which includes the intensity Ii> of the 
j'th template peak. This approach tends to match real 
Qi to the nearest bright peaks in the template pattern, 
naturally accounting for phason shifts, imperfections and 
experimental error that may shift Qi from its ideal value 
and cause an incorrect assignment to some qy with an 
unrealistically small intensity. 

The procedure for finding the best match between the 
real and template patterns, then, is as follows. Choose 
some bright peak in the ideal pattern (No, Mo) with a 
low value of N. An "attempted match" consists of re- 
scaling its magnitude to match the first real peak in the 
pattern, Q\. Identify each remaining peak Q^i in the 
real pattern with the template peak qy which minimizes 
\Qi ~ Next, introduce a goodness-of-fitness pa- 

rameter S\ to characterize the attempted match between 
the {Qi} and the {%'}■ Repeat the process by assign- 
ing (No, Mq) to the second real peak Q 2 and compute its 
goodness-of-fit, S%. After repeating for each real peak, 
use the lowest Sj to decide the best overall match be- 
tween the real pattern and template. Each real peak Qi 
is now assigned the set of indices or, equivalently, 

{N^, M^} of the matching template peak. 

The goodness-of-fit parameter Sj depends on A = 
(\Qi — Qi'\/Qi) where (d) denotes an intensity-weighted 
average (£\ Vli@i)/Cl2i V^i)- A measures the frac- 
tional deviation between Qi and qv and is weighted by 
intensity for the same reasons as above. (Results do not 
change significantly if a/77 is replaced by h or some sim- 
ilar function of intensity.) One challenge is that for each 
good fit of Qi to a template peak labelled by (N, M), 
there is also an equally good fit to the peak r 6k (N, M) 
where k is any integer, due to the self-similarity of the 
lattice; see Eq. (||). This leads to the practical problem 
of deciding which match to choose and also the annoy- 
ance that, unlike a crystal, the indexing is not unique. To 
produce a unique and sensibly standardized indexing, we 
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simultaneously minimize N — (Ni). The goodness-of-fit 
S can be taken to be a linear combination of these two 
parameters, S = aA + N. The results are relatively in- 
sensitive to the choice of a (up to an order of magnitude) 
provided that the A and N contributions to S are both 
non-negligible for known quasicrystals; typical values are 
N w 30 and A rj 0.3, and we used a = 500 in all data 
runs. 

If the real material is a crystal, the best attempted 
match to the quasicrystal template (corresponding to the 
lowest Sj) is still a poor match relative to a true qua- 
sicrystal, so we next introduce statistical tests to measure 
the quality of the match. These tests were found empiri- 
cally by applying them first to known quasicrystals, and 
they involve calculating several quantities for each pat- 
tern. The first two statistical quantities, A and Q± (the 
intensity- weighted average of Q±), are discussed above. 
Quasicrystals have low values of both, representing closer 
matches to brighter peaks. However, finite resolution 
limits the number of peaks present above the noise level; 
while these two parameters clearly separate out FCI qua- 
sicrystals (see Fig. 0), they fail for the SI case, where 1700 
patterns score better than known SI quasicrystals when 
ranking only by these two statistics. We therefore con- 
sider another quantity that involves interrelations among 
the bright peaks in a quasicrystal powder pattern. 

Each peak in a powder pattern has a parity determined 
by r\ = J2i n i- Even (odd) rj corresponds to even (odd) 
parity. Over a finite interval, certain sequences of even 
and odd peaks appear. For a given bright peak at Qq, 
there are other bright peaks at Q = Qo + AQ, where 
AQ = AN + AMt = (AN, AM). All these peaks should 
have low values of Q±, and will therefore be separated 
by small AQ\ = ANt - AM. AM /AN should ap- 
proximate r, the golden ratio, and as is well known, the 
best approximant is a ratio of subsequent Fibonacci num- 
bers. Hence, for every peak in the powder pattern, we 
look for other peaks separated by "Fibonacci intervals," 
(AN, AM), where AN and AM are either successive Fi- 
bonacci integers or constant multiples of these integers. 
In the SI case, we search for two types of intervals: those 
occurring between peaks of opposite parity, which involve 
successive Fibonacci numbers, such as (2, 3) and (34, 55); 
and those occurring between peaks of the same parity, 
which involve four times successive Fibonacci numbers, 
such as (8, 12) and (20,32). In the FCI and BCI cases, 
only even-parity peaks are present, and we seek only se- 
quences with four times successive Fibonacci integers. 

For each ith peak in the real pattern, we count the 
number of other real peaks that are separated from it by 
one of the Fibonacci intervals. If there is only one such 
peak, then the ith peak is considered part of a pair (or 
2-plet); if two, part of a 3-plet; etc. We define M| 
to be 1 if the ith peak is part of an m-plet and zero 
otherwise. Note that M-"^ — 1 implies M| m_1 ^ = 1. For 
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FIG. 1: A plot of the distribution of quantities A and Qj_ for 
eleven patterns identified as FCI quasicrystals in the ICDD- 
PDF (white circles) and 60,000 patterns identified as crystals 
(grey dots). The large square shows the scatter plot in the 
two-dimensional parameter space, in arbitrary units with the 
origin at lower left. Rectangles along either axis show the 
individual histograms of A and Qj_. Means and standard 
deviations for each individual quantity are marked by dashed 
and dotted lines, respectively. The quasicrystals are projected 
to show that they lie at the tail of the distributions. 

example, if the ninth peak in a pattern is separated from 

(2) 

three other peaks by Fibonacci intervals, then Mg = 
M 9 (3) = M 9 (4) = 1 and M<f } =0. The values of M/ m) for 
each peak are combined into intensity-weighted averages 

M (m) = ( M f m) ). 

The three average quantities A, Q± and M*- 4 - 1 are uni- 
fied into a single x 2 statistic, which quantifies the degree 
of quasicrystallinity of a powder X-ray diffraction pat- 
tern. First, the distribution for each quantity, calculated 
for all sixty thousand entries in the ICDD-PDF with a 
dozen peaks or more, is fit to a standard gaussian by 
matching percentiles. Then, \ 2 is calculated from the 
three normalized measures. All patterns with any worse- 
than-average quantity are discarded and the remaining 
patterns with a high value of x 2 represent patterns with 
characteristics most like those of the known quasicrystals. 

Results. The x 2 statistic has proven to be a reliable 
method of identifying and indexing quasicrystal diffrac- 
tion patterns. For example, Fig. [I] shows a scatter plot 
of A vs. Qx for 60,000 patterns from the ICDD-PDF. 
The eleven patterns known to be synthetic quasicrys- 
tals of the FCI structure (by checking the references 
given in the ICDD-PDF) are represented by white circles; 
the remaining patterns are shown in grey. The cluster- 
ing of the eleven quasicrystals apart from the remaining 
points demonstrates the success of the tests in separat- 
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TABLE I: Top quasicrystal candidates for SI, FCI, or BCI lat- 
tices. Type refers to whether candidates were labeled in the 
ICDD-PDF as quasicrystal, synthetic crystal or mineral crys- 
tal. Shown are a typical quasicrystal, the top five synthetic 
and top mineral candidates. The best known quasicrystals 
have x 2 values up to 154.8 (SI) and 397.8 (FCI). 



PDF Num 


Formula 


Lattice 


Type 


x 2 


New Datag 


CdYb 


SI 


QC 


34.7 


27-901 


SnTe 3 8 


SI 


Syn 


28.4 


42-842 


InP 3 


SI 


Syn 


27.9 


44-583 


CaUC-4 


SI 


Syn 


27.4 


21-117 


Cd(Mn0 4 ) 2 -6H 2 


SI 


Syn 


21.4 


38-923 


K 2 NaPdF 6 


SI 


Syn 


19.8 


25-298 


Aktashite, CueHg3As4Si 2 


SI 


Min 


17.9 


48-1437 


Al 68 . 5 Pd 2 2.iMng.4 


FCI 


QC 


92.0 


a a i nc 
4U-100 


Ba3La4o V12U93 


t Ul 


byn 


4U.3 


19-261 


CaYb 2 4 


FCI 


Syn 


38.9 


31-1420 


U0 3 


FCI 


Syn 


34.1 


41-979 


Cro.9Ta5.1S 


FCI 


Syn 


33.8 


27-863 


Sr 7 Yi304(P0 4 )3(Si0 4 )9 


FCI 


Syn 


33.7 


2-691 


Tantalite, (Fe,Mn)Ta 2 6 


FCI 


Min 


25.2 


21-379 


InPtU 


BCI 


Syn 


40.1 


50-1135 


Co 4 Sni 3 Tb 3 


BCI 


Syn 


36.2 


42-1163 


PbioAl 2 F 26 Cl 


BCI 


Syn 


26.8 


45-1164 


AlaoMoLeseTh 


BCI 


Syn 


25.9 


36-1231 


LaNi0 2 


BCI 


Syn 


24.9 


44-1412 


Gratonite, Pb9As 4 Sis 


BCI 


Min 


14.3 



ing out quasicrystal patterns. The few grey dots around 
the white quasicrystal circles are patterns identified as 
crystalline in the ICDD-PDF that might be misidenti- 
fied quasicrystals. However, many of these possibilities 
are eliminated when the remaining test is applied. 

Because no one or two tests is completely effective, we 
combine three quantities into a x 2 statistic. 

Tabic U lists the top quasicrystal candidates for the SI, 
FCI and BCI structures. For each structure, the list in- 
cludes a typical example of a known quasicrystal (except 
BCI, where none is known), the top five synthetic mate- 
rials identified as crystalline, and the top mineral candi- 
date in the ICDD-PDF. The ranking should not be taken 
as an absolute measure. Different versions of the statis- 
tical tests can alter the ranks of some promising candi- 
dates, sometimes significantly, but those shown remain 
among the top-ranked throughout. Our quasicrystal ex- 
ample in the SI case is the recently discovered binary 
alloy, CdYb. || Interestingly, a stoichiometrically similar 
phase, (ostensibly crystalline) Cd 6 Yb, placed seventh in 
the BCI top candidate list. For the BCI case, there is no 
known quasicrystal, and though the top candidates have 
high values of x 2 , the \ 2 values cannot be compared di- 
rectly between different structures since statistical tests 
are different. Nevertheless, the high scores suggest fol- 
lowup studies of the best candidates, particularly InPtU 
and Co4Sni3Tb3. 

A further success is that our automated indexing of 



peaks for the known quasicrystals matches the available 
published indices. In some cases, while the published in- 
dices are aided by two-dimensional diffraction data, our 
automated procedure uses powder data only. Given the 
unique challenges in quasicrystal indexing, described ear- 
lier, the successful indexing suggests that our matching 
of template and real patterns automatically incorporates 
standards obtained by individual analysis. Hence, as 
more quasicrystals are discovered and added to catalogs, 
our procedure could provide a standardized indexing. 

Future work will proceed in several directions. A num- 
ber of materials merit further study, and, intriguingly, the 
most promising examples are in the BCI class where no 
quasicrystal has yet been found. Systematic studies of 
synthetic materials with nearby stoichiometries may also 
be merited. Finally, we would like to collect the powder 
diffraction patterns of other materials, including known 
quasicrystals and crystal approximant phases, to further 
refine our tests. We are interested in collaborating in ex- 
ploring the leading candidates, only some of which have 
been given in Table §. Those interested are encouraged 
to contact PJL and PJS. 

We thank the ICDD for providing a copy of the ICDD- 
PDF; P. Tandy (Natural History Museum, London), J. 
Post (Smithsonian) and L. Cabri (CANMET) for lend- 
ing us mineral samples for pilot studies; and A. -P. Tsai 
for sending us the CdYb powder pattern. PJL thanks B. 
Fujito and M. Yogo for assistance with coding and statis- 
tics, respectively. This work was supported in part by the 
US Department of Energy grant DE-FG02-91ER40671 
(PJS) and NSF Materials Research and Engineering Cen- 
ter grant DMR-94-00362 (NY). 
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